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1. Introduction 


Oscillations between active neutrinos v a and sterile neutrinos v a have been discussed in 
various astrophysical and cosmological environments. Such examples include the oscillation 
between i/ e and v a as solutions to the solar neutrino problem [1], and oscillations between 
v a and u s m Big Bang Nucleosynthesis which may affect the primordial 4 He abundance 
[2-6]. Constraints on neutrino mixing parameters have been obtained by requiring the 
predicted 4 He abundance under these neutrino mixings not to exceed the observed value 
[3-5]. Sterile neutrinos are also suggested as dark matter candidates if they have proper 
mixings with active neutrinos [6]. In this paper we attempt to consider v t -v s mixing in 

another astrophysical environment-the Type II supernova (hereafter we omit “type II” for 
simplicity) core. 

Since there isn’t a fully convincing numerical supernova model at present, we try to 
construct our arguments as model-insensitively as possible. Supernova explosions occur 
when massive stars (> 8M 0 ) reach the end of their evolution [7]. For such massive stars, 
an iron core is developed at the center of the progenitor star, surrounded sequentially by 
shells of silicon, oxygen, nitrogen, carbon, helium, and finally an envelope of hydrogen. 
Nuclear reactions occur at the interfaces of these shells and between the shell and the 
core. The initial iron core of the progenitor star has a density of about 10 9 g/cm 3 , and 
is primarily supported by the pressure of the degenerate electron gas. As the nuclear 
burning of the silicon shell continues, more iron is deposited into the core, until it reaches 
the Chandrasekhar mass (» 1.4M©), when the core cannot be sufficiently supported by 
the electron degeneracy pressure and begins to collapse. In the meantime, due to the high 
density and the high temperature of the core, the deleptonization process c - +p - „ + ^ 
occurs rapidly. When the density reaches 10 11 g/cm 3 , the mean free path of neutrinos l 
(« 10 3 / p 14 cm for neutrinos of ~ 20 MeV, where p 14 is the core density divided by 10 14 
g/cm 3 ) is shorter than the radius of the core i? core - 10- 100km. When the density is higher 
than 10 12 g/cm 3 , the diffusion timescale of neutrinos becomes longer than the timescale 
of the collapse of the core. Neutrinos are thereafter trapped inside the core instead of 
escaping freely. The collapse of the core finally stalls and bounces back when the center of 
the core reaches approximately nuclear density fewxlO 14 g/cm 3 , at which the pressure of 
the degenerate neutron gas stops the gravitational collapse. The bounce creates a shock 
wave in the midst of the core that propagates outwards through the infalling matter. It 
is generally believed that if the core is small (< 1.4M 0 ), the shock wave will have enough 
energy to blow the matter outside the core away and lead to an explosion event [8]. If the 
core is sufficiently massive, however, the shock wave quickly stalls within the core by losing 
energy m dissociating large quantities of infalling nuclei into free nucleons. It is currently 
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believed that in this case neutrinos deposit their energies into the shock several tenths of 

a second after the bounce and reheat the shock wave to eventually blow the outer layer 
away [9J. 

It is estimated that about 99% of the energy is emitted in the form of neutrino bursts. 
Therefore, the role of neutrinos in the supernova event is essential [9], At the end of the 
coUapse, electron neutrinos axe trapped instead of streaming out of the core. The 
and anti-neutrmos are also produced thermally inside the hot core with an average energy 
of several tens of MeV. This leads to a neutrino sphere with an average energy of ~ 20 
MeV trapped inside a dense core of nucleons and heavy nuclei. The number density of 
" e 1S highly su PPr e ssed at the collapsing stage of the supernova due to the existence of 
the degenerate v e gas. Its number density builds up only after significant amount of u c 
escapes from the core after deleptonization. It was shown by Wilson et a 1. that the * and 

'' annihilating at the ne “‘ r ino sphere can deposit their energies in the shock and lift the 
shock front outwards, thus leading to an explosion [10]. Diffusive neutrinos further cool 
down the core (a proto-neutron star) to a neutron star after the explosion, on a timescale 
of about 10 seconds. This diffusive * e flux has been observed by Kamioka [11] and IMB 
[12] in conjunction with SN1987A, with a typical energy of 15-20 MeV and a timescale 

of about 10 seconds. Therefore, the current picture of a supernova event has received 
significant observational support. 

2. Neutrino Mixing in the Supernova Core 

The above picture assumes the standard electroweak theory. If neutrino mixing is 
present, the picture may change substantially. Many literatures have been focusing on the 
constraints on properties of neutrinos [13]. We will focus on the impact of a v e - v , mixing 
on the supernova event, which hasn’t been explored in full detail before [14]. In particular, 
instead of discussing the general cases of the mixing, we will discuss only the case of a 
resonant oscillation between * e and u„ which are relatively less model-dependent. 

In the flavor basis the effective Hamiltonian of the *„ (active)-*, (sterile) mixing system 
m a medium is given by [15] 

H — — L_ ( ~ ^ m2 cos 20 -f- 4 EV Va Sm 2 sin 20 A 

4£ V <§ra 2 sin 20 Sm 2 cos 20 ) ’ W 

where E is the neutrino energy and Sm 2 = m 2 2 - to 2 is the difference of the squares of the 
two vacuum mass eigenvalues xn\ and m 2 . We assume Sm 2 > 0 when the sterile neutrino 
is heavier. Furthermore, 8 is the vaccum mixing angle and V Va is the effective potential of 
the active neutrino v a in the medium. In a stellar environment [16], 
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( 2 ) 


V Ue - V2G F (N e - 0.5 N n + 2 N Vt -f N U(i + N„ r ), 

V Um = V2G F (-0.5N n + N Ve + 2N U " + N Vr ) 

where N t , N„ and N„ a are the number densities of electrons, neutrons, and neutrinos minus 

t e number densities of their anti-particles, respectively. V„ r is similar to V„ except that 

v, and are interchanged in the equation. The effective potential of a sterile neutrino is 

always zero. Based on CP symmetry = -V „, . In the special situation of a supernova 

core, Av„ and N„ r are negligible, since u„ and u r are only produced in vu pairs via neutral 
current processes. 

A resonance between the two mixing neutrinos occurs when the diagonal elements of 

the effective Hamiltonian eq. (1) become equal, i.e„ when the neutrino energy equals the 
resonance energy E r given by 


E r = 


Sm 2 cos 26 
2Vu 


(3) 


ns, e a supernova core, where spherical symmetry is a reliable approximation, the region 
where such a resonance occurs for neutrinos with an energy E r is a spherical shell with 
radius RUE r ). It can be shown that if the sterile neutrinos resonantly converted from 
the active neutrinos escape freely from the core (we will discuss the condition for this 
ater), the usual continuity equation relating the local rate of change of the active neutrino 
number density n„ a to the divergence of the corresponding local flux is modified to 


n, 


V • - 


d<j>?, nt • dn OVLt 

VE,-^-(E r) + E,^(E r) 


dE 


dE 


Pv. 


(4) 


Here, $ and are the partial flux and number density which only account for v a ’s 
going outward with respect to the surface of constant resonance energy E r (outward radial 
irection m our case). Furthermore, P Q is the probability for an active neutrino going 
t rough the resonance to be converted to the sterile neutrino We have neglected terms 
corresponding to v, -» v a convertions and i/ a -> i/ 3 convertions of inward going v a since 
we are only interested in the global „„ loss by integrating eq. (4) over the whole SN core 
to which the sum of the neglected terms does not contribute. 

In general, P a depends on the rate of decoherence of the neutrino flavor density matrix 
caused by non-forward scattering off the background [14] as well as on the adiabaticity of 
the resonance. In the limiting case where the width of the resonance region is much shorter 
t an the mean free path of neutrinos, which is the case we will do our calculations, 


P a (E, 


exp 


f — 7r<5m 2 sin 2 26 


4 cos 26 
4 


V£ r + 


Er 


-u 


(5) 



This is correct to the first order in the change of E r in the resonance region [15]. P a ~ 1 

corresponds to a highly adiabatic resonance whereas P a < 1 represents a highly non- 

adiabatic resonance. On the other hand, when the width of the resonance region is much 

larger than the neutrino mean free path, P a ~ 0.5 [14], which means the resonance is 

always “adiabatic”. Therefore, it can be seen from our following calculations that our 

constraints that exclude sufficiently adiabatic resonant mixings will equally apply to the 

case with large resonance width, as long as the stream freely outside of the resonance 
region. 

To get an order of magnitude estimate of the neutrino loss rate, we can further simplify 

eq. (4) by setting V E r ■ (d$ onl j dE){E r ) fa ( dE r /dr)ncf(E r )E?/E 3 and (dn oal / dE)(E r ) fa 

nf(E r )E r /E (we drop the subscript v a from now on). Here, c is the speed of light, f(E r ) 

is the occupation number at a state with E r , and E is the average neutrino energy which 

for the partially degenerate u e is of the order of its chemical potential p and for other 

neutrino species (and all antineutrinos) is 3 times their temperatures in the supernova 

core. We can then integrate the second term on the r.h.s. of eq. (4) over the SN core to 

get the additional global v a loss rate N Ies and the additional global energy loss rate Q res 
due to the 


due to the resonant conversion, 

'Rcore 

N 


. 2 dE r 
47 rr n 


dr (c + Vr) f(E r )§P(Er)dr, 
E 


Qres ~ 


J f Rcore 9 dE / p \ 3 

, iTr n -jf( c + Vr)f(E r )^ P(E r )dr , 


( 6 ) 

where v T = E r /\VE r \, the speed of the resonance shell for neutrinos with energy E r . 
The integrals on the r.frs. are dominated by integrands at E r ~ E if the resonance for 
neutrinos with energy E ever occurs inside the core. When E r ~ E, we always expect 
\E\«c- dE/dr since otherwise the neutrinos will build up inside or dissipate from the 
core m H^/c ~ 10 4 s, which contradicts our knowledge of supemovae. Therefore v r and 
E are negleted in the following discussion. If we denote iV tot and Q tot as the total number 
and the total energy of u a inside the core, we can approximate eq. (6) by 

V. A CP (Z)\E r =0„ 

Rear, ^ Q ™ RZr~ Qt °'- (7) 

One crucial quantity in the problem is the scale height of E r , \d\n E T / dr\~\ A rough 
estimate for it is ^ 


dhiE r 

-l 

d In N e ff 

-l 

din n i -1 1 Neff 

dr 


dr 


dr f 1 n 


( 8 ) 


, , — — ixv.igjij.o yj x neutrino aensi 

and N eff = V Va /y/2Gp. Without going into detailed models, 

| din n i -1 

> /. 


dr 


(9) 



In the following sections we apply these results to estimate the loss rates in the col- 
lapsing and the cooling phases of supernovae. 

3. Constraints from the Collapsing Phase 

During the collapsing phase of supernovae, v e will encounter a resonance when the 
v e -v 3 mixing satisfies 

6m 2 cos28 = 2v^G F i? r AW = 2t/2 G F E r (N e - 0.5 N n + 2 JV„.). (i 0 ) 

Taking E ~ lOMeV, and the density of the core p ~ 10 12 g/cm 3 to several times 10 1 * g/cm 3 
we find 10 eV < (Sm < 10 8 eV 2 . At this stage the i/ e ’s are already trapped in the SN 
core. The relevant time scale here is the free fall time scale of about 10- 3 s~ ft„„/f0 lcl 

en 1 re, > (0.1c/fi co „)iV tot , i/ e will never get trapped after the central density of the 
core reaches and exceeds the resonant value. The process is illustrated in Figure 1. It 
has been shown that no explosion occurs when the trapped lepton number is less than 0.9 
times the lepton number before the collapse [17]. Therefore we require 

P ( E )\b=o £ 01 for 10 6 eV 2 < Sm 2 < 10 8 eV 2 (H) 

so that no serious leakage of ^ will occur. Because of eq. (5) this corresponds to 


6m 2 sin 2 20 1 d In E r i -i 


* 1 *-<*) SOI- 

(12) 


Since JVeff ~ n , \dlnE r /dr(E)\ 1 > / ~ 10 3 cm from eqs. (8) and (9). Taking E ~ lOMeV, 
we conclude that v t -v a mixings with 

Sm 2 sin 2 20 > 10 -2 eV 2 and 10 6 eV 2 < Sm 2 < 10 8 eV 2 ( 13 ) 

are ruled out so that the supernova core can trap enough lepton number to explode. 

4. Constraints from the Cooling Phase 

For 6m <10 eV , interesting constraints on the v e -v 9 mixing come from the con- 
sideration of the resonance that occurs after the bounce. After the bounce, N e „ = 
N e - 0.5 jV„ + 2N„ t drops below zero very quickly at the outer part of the core [17,18], (In 
general N e = N p ss N n before the collapse, where N p is the number density of protons. 
After about 1/3 of protons are deleptonized into neutrons, 7V eff will be negative.) However, 
at the center of the core, there is a significant build-up of i/ e , N e „ remains positive. Figure 
2 shows roughly the relation between 7V eff and the position inside the core at about 0.5 
seconds after the bounce, as inferred from figure lOf of ref. 18 (similar relations exist for 
prompt-explosion models [17] ). Thus, there exists a central part of the core with a positive 
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V ’' and an ° Uter part of the core with a negative V„ e . Furthermore, outside the core where 
no significant deleptonization has taken place, V„ e is still positive. Therefore two resonance 

regions may form, (for both v t -v, mixings and v t -v, mixings), one inside the core and one 
outside the core, with 


-Neff ~ 0, 


(14) 


if 6m is much smaller than G F EN eff and can be neglected. If we take E ~ 10 MeV at 
the edge of the neutnno sphere, we can see from figure 2 that when 6m 2 < 10 6 eV 2 , both 
Ve and can be ^onantly converted into sterile neutrinos inside the core. 

In this section we consider the resonance region inside the core. In this case, the scale 
height of E r is much smaller than in the previous section due to eqs. (8) and (14). By 
using eq. (10) and n <Hf / 67r 2 , we find for E r ~ ~E, 

3 ir 2 6m 2 


dlnEt 


dr 


~\e) > l 


iV eff 


n 


>/■ 


V2G f E* 


(15) 


The concern is that if the resonance occurs with a sufficient adiabaticity, the intensive 
production of the freely streaming sterile neutrinos will cool down the core too quickly 
compared with the observed cooling time of about 10 seconds. This translates into 

c ^\e v =q 


rtcore {E) > 1 S6C ’’ or equivalently P|^ =0 (£) > 10' 3 

for Rcote ~ 10 to 100km. By using an approximation 

Dl 6m 2 sin 2 29 , _ , . 

~ cos 26 |V£?r| for P « 

we conclude after combining eqs. (15), (16) and (17) that if 

(8m 2 sin 2 2 6\ / Zir 2 8m 2 l \ 


V 


cos 26 


-){ 


0 


>10 


-3 


(16) 


(17) 


(18) 


y/2G F E 

the proto-neutron star will cool down much faster than the observed timescale by losing 
neutrinos. Taking E ~ 50 MeV, l ~ 10 3 cm, eq. (18) excludes v c -v a mixings with 


8m 2 sin 28 > 30eV 2 and Sm 2 < 10 6 eV 2 . (19) 

5. Possible Tightening of the Bound by Reconversion 

Since the effective potential K. goes from negative at the edge of the core to positive 
outside the core, at least the P e -P s will encounter a second resonance if 6m 2 < 10 6 eV 2 . 
The density outside the core is quite low compared with that of the core. As a result, 
this outer resonance could be highly adiabatic since |rfln E r /dr\~' could be much larger.’ 
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Then all the P, converted from P. deep in the core will be completely reconverted into P. 
which then freely escape from the supernova because the outer resonance is located above 
the neutrino sphere. But these neutrinos have energies >50 MeV as they are originated 
from P.’s deep inside the neutrino sphere. Since the cross section of neutrinos in the u-e 
scattering experiments (Kamiokande and IMB) scales roughly with E 2 , where E is the 
energy of the neutrinos, the fact that we didn’t see any p. with E > 50 MeV in these 
experiments will put a tighter bound on mixings than that of the previous section. 
The average energy of v t observed by Kamiokande and IMB is ~ 20 MeV [13], Therefore, 
the above consideration could extend eq. (16) by two order of magnitude, because even if 
the energy loss rate due to the resonant conversion is only one tenth of the standard loss 
rate, ~ 5 of the 20 neutrinos seen by IMB and Kamiokande should have had an energy 
of >50 MeV. More generally, if in the future a neutrino detector detects M supernova 

neutrinos with none of them having an energy > 50 MeV, the excluded parameter region 
would extend to 

8m 2 sin 20 > ^ eV 2 and 8m 2 < 10 6 eV 2 . (20) 

With the current data, M = 20, the excluded region is 

8m 2 sin 20 > leV 2 and 8m 2 < 10 6 eV 2 . (21) 

One has to keep in mind, however, that eqs. (20) and (21) are only true if the outer 
resonance is highly adiabatic, which for eq. (21) requires the density scale height at the 
outer resonance to be larger than ~ (10/sin2<?)m, which seems fairly reasonable when 
compared with the core dimensions. 

6. Discussion and Summary 

As we mentioned in the previous section, the above calculations are valid only if the 
sterile neutrinos have a mean free path larger than the radius of the core after passing 
the resonance. Outside the resonance region, the mean free path of sterile neutrinos is 
roughly 1/ sin 2 0 m [3-5,14], where 6 m is a typical medium mixing angle. In the case of the 
constraints from the collapsing phase discussed in section 3 this is of the same order as the 
vacuum mixing angle 9 whereas in the case of section 4 due to medium effects this is even 
suppressed compared to 0. Our calculations are thus valid as long as 

sin 2 9 < //tfcore ~ 10“ 3 . (22) 

Figure 3 shows the excluded region on the v t -v s mixing parameter space from our 
consideration of a resonant v e -v 3 mixing in both the pre-bounce and the post-bounce 
supernova core, eqs. (13) and (19), as well as the possible tighter bound eq. (21), and 
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eq. (22) under which the above bounds are valid. It also shows the parameter space 
excluded from Big Bang Nucleosynthesis, which gives a stricter constraints for 6m 2 < 10 3 
eV [3 5]. Tha hatched region in Fig. 3 shows the required mixing parameters for the 
sterile neutrino to contribute between 30% and 100% of the critical density of the universe 
[5,6] (assuming a Hubble constant between 50 km/sec/Mpc and 100 km/sec/Mpc). A 
significant part of the hatched region (at which the sterile neutrino serves as warm dark 
matter) is excluded by our constraints from the supernova. 

In summary, we have considered the effect of a resonant v t -v B mixing on a type 
II supernova. By requiring that a supernova retains enough leptons to explode and is 
not cooled down within a fraction of a second, we obtain constraints on the u e -i/ 9 mix- 
ing parameter space. The excluded region includes the region required for the sterile 
neutrino to be the warm dark matter through v t -v a oscillations in the early universe. 
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Figure Caption: 

Figure 1. An illustration of the configuration of V„, inside the core during the core collapse 
With (the dashed line) and without (the solid line) the i mixing that satisfies eq. ( 13 ). 

Figure 2. An illustration of N e g inside the core at the bounce and several tenths of a 
second after the bounce. R n , is the radius of the neutrino sphere. 

Figure 3 . On the u.-u. mixing paramter space, the area enclosed by the solid line is 
excluded by our supernova consideration. The area enclosed by the long-dashed line is 
the possible extended bound due to the resonance outside the core. The area enclosed by 
short-dashed line is excluded by Big Bang Nucleosynthesis. The hatched region shows the 
parameters for a sterile neutrino dark matter candidate (generated a v c -v. mixing in the 
early universe) to contribute 30 % to 100 % the critical density of the universe today (the 
Hubble constant is taken to be between 50-100km/sec/Mpc). 
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